Dynamical simulation of transport in one-dimensional quantum wires 
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Transport of single-channel spinless interacting fermions 
(Luttinger liquid) through a barrier has been studied by nu- 
merically exact quantum Monte Carlo methods. A novel 
stochastic integration over the real-time paths allows for di- 
rect computation of nonequilibrium conductance and noise 
properties. We have examined the low-temperature scaling of 
the conductance in the crossover region between a very weak 
and an almost insulating barrier. 

PACS numbers: 72.10.-d, 73.40. Gk 



Studies of transport in one-dimensional (ID) interact- 
ing Fermi systems have gathered novel physical insights 
and led to the development of several new techniques 
over the past few years. Instead of the usual Fermi liq- 
uid theory, the low-temperature properties of an ideal- 
ized ID quantum wire are described by the Luttinger 
liquid model [|l]j2). Remarkable transport behaviors are 
expected for a Luttinger liquid in the presence of impu- 
rities or barriers, and several interesting theoretical 
and experimental studies have emerged recently. In 
this Letter, we describe a real-time quantum Monte Carlo 
(QMC) method that enables us to directly address the 
dynamics of a Luttinger liquid, and use it to investigate 
the low-temperature transport properties of a ID quan- 
tum wire. 

So far the only numerical approach to this problem has 
been the one given by Moon et al. |Q . These workers ex- 
amined transport in a Luttinger liquid using Euclidean- 
time QMC simulations by analytically continuing numer- 
ical data to real time with Fade approximants. Math- 
ematically, this analytic continuation is ill-posed and 
small statistical errors in the imaginary-time data could 
be magnified immensely |l0[| . Physically, imaginary- 
time data collected at low temperatures contain pre- 
dominantly ground-state information whereas real-time 
dynamics is controlled by the excitation spectrum. Al- 
though the findings of Ref. |^ are certainly reasonable, 
a direct real-time simulation which avoids the possibly 
troublesome analytic continuation is desirable. Further- 
more, a real-time simulation needs not rely on the Kubo 
formula and allows for direct computations of nonequi- 
librium transport and noise properties beyond the reach 
of previous methods. 

Common to all real-time QMC simulations is the ubiq- 
uitous and fundamental dynamical sign problem [p^-p^. 
It arises in the stochastic summation over the system 



paths when the real-time propagators are oscillatory. 
The quantum-mechanical interference between different 
paths leads to a vanishing signal-to-noise ratio at very 
long real times and the simulation becomes effectively 
unstable. However, it is generally possible to treat 
intermediate-to-long times by employing a partial sum- 
mation scheme p4[ which reduces the effective number of 
variables subject to the Monte Carlo sampling. If guided 
by physical intuition, such a partial summation scheme 
can largely circumvent the dynamical sign problem. Pre- 
vious applications to dissipative tight-binding models in 
the quantum-chemical context 
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15 1 have demonstrated 
the practical usefulness of this concept. Here, we adapt 
this idea to transport in interacting ID quantum wires, 
and as we demonstrate below, the method is able to 
extract numerically exact results in the interesting low- 
temperature region. 

Transport in Luttinger liquids depends crucially on 
a dimensionless parameter g which characterizes the 
electron-electron interaction strength. For the non- 
interacting case, g = 1. For g < I, interactions are 
repulsive, and a barrier of any height leads to a perfect in- 
sulator at zero temperature. At asymptotically low tem- 
peratures, Kane and Fisher (KF) ||] have predicted the 
following scaling behavior for the linear dc conductance: 



G - T^/f- 



(1) 



Most analytical methods applied so far assume ei- 
ther a very small or a very large barrier height (impu- 
rity level) Vq compared to the fermion bandwidth ujc- 
For instance, perturbative arguments lead to Eq.(^ in 
the almost insulating case Vq/oJc ^ 1, while exact re- 
sults find the same scaling law for the opposite 
limit Vq/lOc ^ 1. Unfortunately, little is known in the 
crossover region between the two. A recent study |l6| has 
shown that the usual bosonized impurity Hamiltonian 
can be used to describe this crossover and is valid even 
for ttVo/cJc > 1, which corresponds to an impurity level 
outside the fermion band. Interestingly, Guinea et al. 
|]l7| have suggested a different low-temperature scaling 
form for the weak-barrier case 



(2) 



While exact results give ci = for both Vq /uJc ^ 

1 and Vq/oJc ^ 1, they cannot rule out the scaling 
term for intermediate barrier heights ttVo/lUc ~ 1. A 
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naive reasoning leading to Eq. (|2|) for intermediate barri- 
ers is as follows. While the anomalous KF scaling expo- 
nents hold near the quantum critical points Vq/wc — 
and oo, scale invariance may not survive far away. One 
might then expect regular analytical behavior G ^ T" 
with even power n, and the lowest-order term would give 
as shown in (^. If this scaling form is indeed true, at 
sufficiently low temperatures the scaling would dom- 
inate for g < 1/2 and the KF law (|^) could only be 
observed for g > 1/2. This prediction is clearly of practi- 
cal importance for future experiments on quantum wires. 
As an application of our QMC method, we will address 
the validity of (||) here. 

The action for the Luttinger liquid transport problem 
is usually formulated in terms of a phase field 6* Q . In the 
presence of a short-ranged impurity potential of effective 
strength (barrier height) Vq, the (real-time) Lagrangian 
density takes the form j|] 

C^{l/2g)d^9{x,t)d''e{x,t) (3) 
+Vq d{x) cos[2\/^6l(x, t)] + (eV/y/^) 5{x) e{x, t) , 

where V denotes an external static voltage. Since one can 
eliminate all degrees of freedom away from the barrier by 
Gaussian integrations, a naive real-time QMC simulation 
scheme would have to perform a stochastic summation 
over all possible 9{x — 0, t) paths. After discretization of 
the time line, one is thus left with the task of sampling 
a high-dimensional integral where each of the integra- 
tion variables 0{O,tj) corresponding to time slice j is a 
c- number. 

For numerical reasons, it is now advantageous to per- 
form the real-time simulation within an equivalent dis- 
crete model (a "charge representation"). Discrete vari- 
ables significantly reduce the relevant configuration space 
subject to Monte Carlo sampling compared to the 6 vari- 
ables which are continuous. Their use constitutes the 
first major step of our partial summation scheme. One 
can generally transform Luttinger liquid transport prob- 
lems involving elastic scattering by short-ranged impuri- 
ties into a discrete representation. This is achieved by 
expanding the impurity propagator [p^ , 

dt' COs{2y/^e{0,t')) = 



exp 



E 



dtj^i 



(4) 



^exp 



2V^i I dt' qj{t')e{0,t') 



where qj{t') — J2i=i ~ ^j) is built from the discrete 
variables = ±1. 

The discrete "charges" can be understood as 
Hubbard-Stratonovich fields. Introduction of these aux- 
iliary fields permits us to perform the now Caussian in- 
tegration over the Luttinger liquid bosons 9. Naturally, 



such mapping onto a discrete charge representation is not 
confined to Luttinger-liquid-type electron-electron inter- 
actions and can be done for any Coulomb interaction 
potential. Furthermore, one can use similar charge rep- 
resentations based on (^ for a multi-barrier problem, 
allowing for studies of resonant tunneling or related phe- 
nomena. In the following, we will focus on Eq. (|^) , which 
describes a Luttinger liquid impinging on a single barrier. 

The nonlinear conductance G{T,V) — 81/ dV for 
transport of the Luttinger liquid through a barrier can 
be readily computed from the bosonized current opera- 
tor, J = (e/y/n) 9{0). In fact, the conductance is 



G(T, V) = (e/V^) (d/dV) lim (^(0, t)) 



(5) 



where the average is taken using (^. In practice, the 
expectation value (^(0, t)) reaches a plateau value after 
some time t*, and if the real-time QMC scheme is able 
to reach this plateau, the conductance can be computed 
directly from the plateau value. Since t* grows very large 
at extremely low temperatures, our method becomes in- 
creasingly costly near zero temperature. However, as 
shown below, the dynamical sign problem can be effec- 
tively circumvented to allow for studies in the interesting 
low-temperature scaling region. 

Evaluating (||) by QMC requires construction of a dis- 
cretized path-integral expression for {9(0, t)), where t is 
the maximum length of the time line considered in the 
simulation. Eliminating all 9 degrees of freedom away 
from the barrier, (^(0, t)) is expressed as a weighted 
sum over all sets of forward-backward real-time paths 
{9f{0,t'),9b{0,t')} with t' going from f = to t. In 
the next step, we switch to the charge representation us- 
ing Eq.(^). One can then integrate out the 9(0) paths 
as well at the price of introducing the auxiliary discrete 
charge paths qf{t') and qb{t'). The resulting effective ac- 
tion contains nonlocal couplings between {qf,qb\ due to 
the eliminated 9 modes. 

Next we express qf and qh in terms of their sym- 
metrized and antisymmetrized linear combinations. This 
constitutes the second major step with regard to the par- 
tial summation scheme. The symmetrized component 
corresponds to a quasiclassical degree of freedom since it 
describes propagation along the diagonal of the reduced 
density matrix. The antisymmetrized path, on the other 
hand, describes quantum fluctuations corresponding to 
the off-diagonal elements of the reduced density matrix. 

The basic idea of the partial summation scheme is to 
split off part of the Monte Carlo problem to be evalu- 
ated non-stochastically, while leaving the remainder for 
the stochastic summation. This can be achieved by 
blocking states together in the course of the Metropolis 
walk: instead of sampling states individually, one non- 
stochastically probes the local surroundings of a given 
state to account for its degree of destructive interference 
[njl . Such a blocking strategy avoids regions with highly 
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destructive interference which do not contribute to the 
path integral in the first place but are responsible for 
the dynamical sign problem. In the present case, this 
strategy can be implemented very efficiently by noting 
that one can sum over all symmetrized components ex- 
actly. This is possible since there are no quasiclassical 
self-interactions in the effective action. In the end, only 
the antisymmetric paths are left for the Monte Carlo. 

Any antisymmetric path can now be characterized by 
"kinks" fji and the respective kink times ti. In the QMC 
calculation, we have to discretize time into N sufficiently 
small slices of length 5t — t/N, and no more than one 
kink is allowed on each time slice. The possible kink 
values are yi = 0,±1,±2 corresponding to the charges 
^ = ±1 of the original forward-backward charge paths. 
The Monte Carlo algorithm then performs a stochastic 
summation over all possible kink configurations {y}. We 
have used both kink migration and double kink flip moves 
to generate the Metropolis trajectory. The MC weight 
function is defined such that {9(0, t')) is evaluated only 
at the endpoint t' = t. Our algorithm channels all efforts 
into generating relevant configurations at this one time 
point alone. 

To check our QMC simulation scheme, we have carried 
out two independent tests: (a) short-time exact enumer- 
ation results are accurately reproduced, (b) additional 
QMC simulations for the static noise (/^) show that the 
familiar Johnson-Nyqvist formula is fulfilled at all tem- 
peratures. The results presented below employed dis- 
cretizations on the order of VoSt ~ 0.3, and an upper 
time limit of Vat = 20 was sufficient to reach the plateau 
region of {9{0, t)) for the temperatures studied here. Typ- 
ically, numerical results for the low-temperature region 
were obtained using 25 to 250 million passes and our 
code performs at an average speed of 1 CPU hour per 
miUion passes on an IBM RISC 6000 Model 590. 

Next we present numerical results for transport of a 
Luttinger liquid through a barrier. Figures Q and |^ 
show simulation data for the linear conductance G at 
the two interaction strengths g — 1/S and g — 2/3, re- 
spectively. We have studied two different barrier heights, 
Vq/llJc — 1/6 and Vq/uJc = 1/3. The value g — 1/3 is 
of direct relevance to tunneling of edge state excitations 
in the fractional quantum Hall effect , while recent 

transport experiments in quantum wires |^ have been 
carried out at g = 0.67 ± 0.03. Our results are of imme- 
diate interest for these experiments. 

The main panels show the QMC data over a wide range 
of temperatures. We compare our data to the thermo- 
dynamic Bethe ansatz solution by Fendley, Ludwig, and 
Saleur (FLS) §] which is exact in the weak-barrier case 
Vo/lOc ^ 1. The Kondo temperature is defined as |^ 



with ti 



/3 



1.19599 and t 



2/3 



0.80549. The scafing 



X = dg{TB/Ty-' , 

where dg is readily determined from the golden rule limit 
§ (specifically, di/3 = 0.74313 and ^2/3 = 0.99806). 
The FLS solution reproduces the golden rule behavior 
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FIG. 1. QMC data for the linear conductance G{X) at 
g — 1/3. Barrier heights are Vq/Wc = 1/6 (circles) and 
Vo/iOc = 1/3 (diamonds). The solid curve is the FLS solu- 
tion |5j. Inset: Low-temperature behavior. Error bars denote 
one standard deviation, and the solid curve is again the FLS 
solution. 
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FIG. 2. QMC data for G(X) at g = 2/3. Barrier heights 
are Vb/ti'c = 1/6 (circles) and Vb/ti'c = 1/3 (diamonds). The 
solid curve is the FLS solution, the dashed curve is a guide to 
the eye only. Inset: Low-temperature behavior. 

G/(eV/i) = 5(1 - X^) at T > Tb (small X), while for 
T < Tb it matches the KF scaling law (p. 
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G{X)/{^/h) 



-2/3 



(large X) , 



(6) 



where K1/3 = 3.3546 and K2/3 = 0.8177. 

For g = 1/3, we compare our numerical results with the 
FLS curve at two values of the barrier height, Vq/uJc — 
1/6 and 1/3 (see Fig. |l|). For the smaller barrier, we 
reproduce the Bethe ansatz results well. This is quite re- 
markable since one might have expected more deviations 
from the FLS solution which assumes Vq/wc ^ 1. In the 
FLS solution, the scaling region is found only at very 
low temperatures T/Tb below w 0.15. This is corrobo- 
rated by our results which also indicate the absence of 
asymptotic contributions to the linear conductance. 
The low-temperature KF scaling ([T]) thus seems to hold 
even for the intermediate barrier heights studied here. 
Over the higher temperature range 0.15 < T/T-q < 0.4, 
an apparent scaling is reached which, however, does 
not persist into the asymptotic low-temperature regime. 
In the case of the higher barrier, rather large deviations 
from the FLS solution are seen, especially at low tem- 
peratures. Remarkably, the experimental data in Ref. ||] 
exhibit the same qualitative deviations. 

In Fig. |[ data are presented for g 2/3 correspond- 
ing to the recent quantum wire experiments by Tarucha 
et al. Even for the intermediate barrier heights 

Vq/uJc — 1/6 and 1/3 considered here, we observe the 
G ^ T scaling (^. While our data fully reproduce the 
FLS curve for Vb/^^c = 1/6, there are significant devia- 
tions for a higher barrier. For Vq/ujc — 1/3, we find a 
smaller value ^^3/3 = 0.477 ± 0.005 instead of the FLS 
value 0.8177, but the low-temperature KF scaling (|^) 
still holds. The intermediate-barrier-height effects for 
(7 = 2/3 and 3 = 1/3 are qualitatively similar. 

In conclusion, we have developed and applied real-time 
quantum Monte Carlo simulations to low-temperature 
transport of a Luttinger liquid through an arbitrarily 
high barrier. The close agreement of our data with ther- 
modynamic Bethe ansatz results is rather remarkable 
since we have studied impurity levels near the fermion 
band edge. For Vq/uJc = 1/6, the FLS solution is 
completely reproduced, and only in the higher barrier 
Vq/lOc = 1/3 case are significant deviations observed. 
Recent speculations about possible contributions to 
the linear conductance for intermediate barrier height are 
determined by our calculations to be incorrect, and the 
low-temperature T^/^"^ scaling is supported for the full 
crossover from a weak to an almost insulating barrier. 
We believe that the method presented here will be use- 
ful in future studies of nonequilibrium transport, noise 
properties, and resonant tunneling. 
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